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Collective charge density excitation spectra of both a spatially separated 
two-component quasi-one dimensional (ID) quantum plasma, as existing in 
a semiconductor double quantum wire structure, and a ID homogeneous 
electron-hole plasma, as appropriate for a photoexcited semiconductor quan- 
tum wire system, are calculated within the two-component random-phase- 
approximation. We find two phase fluctuation collective modes, one (optical 
plasmon, OP) with energy proportional to q\ ln(qa)| 1//2 is the total in-phase 
(out-of-phase) charge density oscillation of the system and the other (acoustic 
plasmon, AP) with a linear energy dispersion as q — ► is the neutral out-of- 
£NJ ! phase (in-phase) charge density oscillation of the system for the situation 

where the two components have the same (opposite) charges, where q is the 
00 ■ ID wave vector and a is a characteristic ID confinement size. In contrast to 

higher dimensional systems we find the neutral long wavelength AP mode to 
be generically undamped by Landau damping effects due to the severe sup- 
pression of single particle excitations in ID systems. We also investigate the 
[ effect of impurity scattering on the collective mode dispersion and damping, 

and calculate the collective mode spectral weight by obtaining the dynami- 
cal structure factor. We find that both OP and AP modes are overdamped 
by impurity scattering below some critical wave vector. We find that in the 
• i-H ! long wavelength limit the spectral weight is carried mostly by the OP, but 

^ | the spectral weight of the undamped AP mode at finite (but not too large) 

wave vectors is comparable with that of the OP mode, making it viable to 
observe the AP mode in semiconductor quantum wire systems. The effect of 
the interwire electron tunneling in a biwire system on the collective charge 
density excitation spectra is also studied. We discuss the mode dispersion and 
damping from an effective Luttinger liquid perspective as well and in some 
cases include local field corrections in our calculations. 
PACS Number : 73.20.Mf; 71.35.Ee; 71.45. Gm; 71.45.-d 
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I. INTRODUCTION 



Recent advances in fabrication techniques involving molecular beam epitaxy and lithog- 
raphy have now made it possible to make narrow GaAs-based quasi-one dimensional (ID) 
electronic systems with lateral dimensions of the order of the Bohr radius In these so- 
called quantum wire structures, the motion of charge carriers is confined in two transverse 
directions but is essentially free (in the effective mass sense) in the longitudinal direction. 
Study of collective modes in reduced dimensional electron systems in semiconductor nanos- 
tructures is a subject of growing experimental and theoretical interest. Experimentally, 
far-infrared optical spectroscopy and resonant inelastic light scattering spectroscopy 
|§-|| have been used to study quasi-lD elementary electronic excitations. Several theoret- 
ical studies, mostly based on the random phase approximation (RPA), have been reported 
on the energy dispersion of elementary excitations in semiconductor quantum wires. The 
measured ID intrasubband plasmon dispersion agrees remarkably well with the RPA predic- 
tions 0. The quantitative agreement between RPA predictions and the experimental results 
H was later explained by the fact || that the calculated RPA plasmon dispersion and the 
Tomonaga-Lutinger theory for the collective charge density excitations of the ID electron 
system are equivalent at long wavelengths by virtue of the vanishing of vertex corrections 
|J in the irreducible polarizability of the ID electron gas. 

The two-component quasi-lD systems which are composed of electrons and holes have 
been generated in a wide variety of semiconductor quantum wires by optical pumping. 
Spatially separated two-component quasi-lD systems can be produced by additional lat- 
eral confinement to produce a double quantum wire structure in GaAs-GaaAli_ x As double 
quantum well systems. The collective modes of the two-component plasma also play impor- 
tant roles in the many-body physics of carriers such as screening of the Coulomb interaction 
potential and in Coulomb drag problem. In this paper, we calculate the dispersion and the 
spectral weight of longitudinal collective modes (with and without spatial separation) in 
quasi- ID two-component plasmas as in semiconductor double quantum wire structures or in 
photoexcited homogeneous ID electron- hole plasmas (EHP). Our theory is based mostly on 
the two-component RPA (which should be well- valid in ID systems); in the ID EHP case 
we go beyond RPA to include local field corrections to assess the validity of ID RPA. The 
collective modes of the two-component plasma have been widely studied in higher dimension 
1^-12 1, and recently, numerically in quasi-lD EHP ||13|| . 

It is well-known that a two-component plasma has two branches of longitudinal collective 
excitation spectra called the optical plasmon (OP) and the acoustic plasmon (AP), where 
the density fluctuations in each component oscillate in-phase (OP) and out-of-phase (AP) 
respectively relative to each other, assuming the two charge components to be the same. 
As expected, we find two plasma modes in the two-component ID system; the in-phase 
OP with a typical ID plasmon dispersion with the frequency proportional to q\ \n{qa)\ 1 ^ 2 at 
long wavelengths and the out-of-phase AP with its frequency linear in q as q — > 0, where q 
is the ID wave vector and a is a characteristic ID confinement size (i.e., roughly the wire 
width). Unlike 2D and 3D systems, where the electron- hole pair continua have only an upper 
boundary, qvF + q 2 /2m, the ID pair continuum has also a low energy gap under the lower 
boundary \qvF~ q 2 /2m\, implying no low energy single particle excitations (SPE) are allowed 
in ID. This gap, the non-existence of low energy SPE in ID Fermi systems, arises from 
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severe phase space restrictions (Pauli exclusion principle) imposed by ID energy-momentum 
conservation, and leads to the possibility of the existence of an undamped plasmon mode in 
multi-component ID systems H] . We find a low energy undamped AP mode in the double- 
quantum- wire system for arbitrary interwire separations (including the zero separation case). 
In general, in higher dimensions, the collective mode spectral weight is mostly carried by 
the OP, which makes the experimental observation of the AP mode in higher dimensions 
rather difficult. In contrast, we find that in quasi-lD systems the spectral weight of the 
undamped AP mode is comparable to that of the OP mode and the spectral weight of the 
AP mode increases as the spatial separation increases, which is similar to the corresponding 
2D situation. The greater relative spectral weight of the AP mode in ID is directly related 
to the non-existence of low energy ID SPE. 

The plasma modes in ID systems have vanishing frequencies at long wavelengths, and in 
general their energies at experimentally accessible wave vectors are small (~ 1 meV). Their 
observability, therefore, depends rather crucially on the absence of appreciable damping or 
broadening in the system. While the absence of long wavelength Landau damping (to SPE) 
makes it feasible to discuss the observability of the AP mode, one must also have very small 
impurity broadening in order for the plasmon modes to be observable. Impurity scattering 
causes the carrier motion to be diffusive assuming the scattering to be weak. In principle, any 
impurity disorder, no matter how weak, localizes all single electron ID states; for weak dis- 
order, however, the localization length is larger than the system size and electron dynamics 
could be considered to be diffusive. The diffusive nature of the electronic dynamics strongly 
affects the plasmon dispersion in the long wavelength limit because the plasma frequency 
is vanishingly small at long wavelengths. The ID plasmons therefore become overdamped 
due to impurity scattering induced level broadening at sufficiently long wavelengths. This 
plasmon level broadening may actually play a significant role in determining the physical 
characteristics in quantum wires. The collisional damping of the plasmon modes due to 
impurity scattering occurs in addition to the Landau damping which arises from the decay 
of the plasmon to SPE electron-hole pair excitations. (As emphasized above, while Landau 
damping is suppressed in ID, impurity broadening is always present unless the system is 
perfectly pure.) In the single-component plasma, it is known that, as the impurity scat- 
tering increases, the collisional damping becomes strong and the plasmon mode becomes 
overdamped and disappears at small q in ID and 2D systems (by contrast, this is a small 
effect in 3D since the plasmon energy is finite at zero wave vector) pTSfl . Inclusion of impu- 
rity scattering in ID and 2D overdamps the plasmon below a critical wave vector |TB|.|TB 



This disappearance of the plasmon spectral weight at small q has important consequences 
for the existence and characteristics of the Fermi surface in a one dimensional electron gas 
I6fl . In this paper we investigate the plasmon damping of the two-component ID plasma in 



the presence of impurity scattering. (As noted above any finite impurity disorder leads to 
Anderson localization of ID one electron states, and therefore the ID electron dynamics in 
the presence of impurity disorder is exponentially localized, not diffusive - this is, however, 
a purely academic point of little relevance to high quality GaAs quantum wire systems with 
low disorder, where the localization lengths are many microns long and are in fact longer 
than the ID sample lengths |I7||; we are therefore justified in treating the effective elec- 



tron dynamics as being diffusive, neglecting the Anderson localization effects in high quality 
quantum wires of relevance to plasmon experiments.) 
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We also discuss the effect of interwire electron tunneling on the bi-wire collective mode 
spectra. The most important qualitative feature of the plasmon dispersion in the presence 
of interwire tunneling is that, even though the in-phase OP mode u + depends weakly on 
the tunneling, the out-of-phase AP mode u;_, which is purely acoustic in the absence of 
tunneling (cj_ ~ q), develops a plasmon gap at q = in the presence of nonzero tunneling; 
the plasmon gap in the AP depends nontrivially on the ID electron density n, the interwire 
tunneling amplitude t, and the interwire distance d. Tunneling, therefore, effectively changes 
the AP (uj_ ~ q) into an optical mode with a long wavelength gap. 

We use the RPA for most of our calculations. As emphasized, RPA is a very good 
approximation for the collective mode dispersion in ID quantum wire structures by virtue 
of the essential vanishing of all vertex corrections to the ID polarizability function. We 
consider the zero temperature case and assume that our quasi-lD system has an infinite 
square-well confinement with a finite width a in the y-direction and a zero width in the 
z-direction. It is easy to include a finite width in the z direction but our results do not 
change qualitatively. In typical quantum wires, the wire width along the growth direction 
(z-direction) is much smaller (by an order of magnitude) than that in the lateral direction (y- 
direction), justifying our 2D finite width model for the ID wire. In the finite-width 2D model 
we use a square well confining potential with infinite barriers at y = —a/2 and y = a/2, 
with the electron dynamics along the a;-axis being free. Confinement of the electrons in 
the y and z directions leads to the quantization of energy levels into different ID subbands, 
and we assume the ID extreme quantum limit where only the lowest subband is occupied 
by carriers, and all the ID higher subbands are neglected. We study the plasmon modes 
of the spatially separated two-component system by taking two identical parallel quantum 
wires separated by a distance d. We include results for the case d = because it applies 
directly to the photoexcited ID EHP where electrons and holes are not spatially separated. 
We consider two different spatial configurations for our two-component bi-wire structure: 
two ID wires with a spatial separation d being parallel to each other aligned in the (1) x — z 
plane (separated along the y-axis) and (2) x — y plane (separated along the z-axis). Since the 
difference in our calculated results between these two configurations is within one percent, 
we only give the results for the configuration (1). 

It is well-known that ID electron systems are fundamentally different from their higher 
dimensional counterparts because ID interacting electrons form a Luttinger liquid, not a 
Fermi liquid, and any electron-electron interaction, no matter how weak, destroys the Fermi 
surface, i.e., the discontinuity in the momentum distribution function, in ID. It turns out, 
however, that this distinction, while being profound at a fundamental theoretical level (be- 
cause it implies the non-existence of quasiparticles in ID systems), is rather irrelevant to the 
understanding of the collective mode spectra and their experimental realization in semicon- 
ductor quantum wire structures. As has been emphasized elsewhere, RPA gives a very good 
account of the collective mode properties in ID semiconductor quantum wire structures, and 
we discuss this issue in some details in this paper in the context of the two-component ID 
plasma. 

The rest of the paper is organized as follows: In Sec. II we provide the basic RPA 
formalism as well as the calculated analytic and numerical results for plasmon dispersion, 
damping, and spectral weight by obtaining the dynamical structure factor (which is a direct 
measure of the light scattering intensity); in Sec. Ill we study the effects of the interwire 
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electron tunneling in a bi-wire system on the collective charge density excitation spectra; in 
Sec. IV the plasmon dispersion and damping in the photoexcited ID EHP is studied in both 
RPA and including local field corrections; in Sec. V we briefly discuss the acoustic plasmon 
mode of the two-component ID Luttinger liquid; we conclude in Sec. VI providing a brief 
summary of our results. 



II. RPA THEORY AND RESULTS 

In general, the density fluctuation spectra or the plasmons or the longitudinal collec- 
tive modes of a solid state plasma are given by the poles of the density-density correlation 
function, or equivalently by the zeros of the dynamical dielectric function. For a multi- 
component plasma the collective modes are given by the zeros of the generalized dielectric 
tensor e. In the two-component spatially separated ID plasma the generalized dielectric 
tensor with wire indices i, j (— 1, 2 indicating the two quantum wires) is given (neglecting 
any inter-wire tunneling) within the RPA by 

€ij{q, u) = 5ij - \ u (q)\\ JJ (q. u), (1) 

where q is the ID wave vector in the x-direction of free motion, u is the mode frequency (H = 
1 throughout this paper), Vij{q) the Coulomb interaction in the ID wire representation, and 
U.jj(q, u) the ID non-interacting irreducible polarizability function for the j-th "component" 
(i.e., the j-th wire with j = 1,2 corresponding to the two wires). In an ideal ID electron 
system the Coulomb interaction V(q) in the wave vector space is logarithmically divergent, 
but in our realistic finite width quantum wire model with the confinement width a we 
obtain the intra-wire [inter-wire] Coulomb interaction matrix element V(q) [U(q)} by taking 
the quantizing confinement potential to be of infinite square well type [[T6|j , 

V(q) = V 11 (q)=V 22 (q) = — f(q), (2) 



U(q) = V 12 (q) = V 21 (q) = ±—g(q), (3) 

where e& is the background high-frequency lattice dielectric constant, the + (— ) sign in 
Eq. (|]) refers to interaction between the same (opposite) types of charge in the two wires, 
(the two wires could in principle have electrons and holes respectively, rather than being 
just all electrons or all holes in both wires), and f(q), g{q) are form factors associated 
with confinement matrix elements, which are given, for our infinite square well confinement 
model, by 

f{q) = f dxK (\qax\)h(x), (4) 
Jo 



g(q) = dxK (qyj (ax) 2 + d 2 ^j h(x), 



(5) 



for two wires parallel to each other aligned in the x — z plane. Here K$ is the zeroth-order 
modified Bessel function of the second kind, and the function h(x) = (1 — x) [2 + cos(27nr)] + 
3 sin(27nr)/(27r). If the two wires are parallel to each other aligned in the x — y plane the 
inter-wire Coulomb form factor becomes [with f(q) remaining the same] 

g(q) = / dxK (\q(ax + d)\) h(x). (6) 
Jo 

As q — > 0, f(q),g(q)~K (\qa\) + C, where C is a constant (independent of d) for f(q), but 
a function of the inter- wire spatial separation d for g(q). Note that as d — > U(q) becomes 
V(q), as it should, because the two wires are identical. 

The irreducible ID polarizability Tlo(q, to) in the pure ID system is given by the bare 
bubble diagram. The analytic form of 11° (q, uj) for complex frequency is given by [/J 

uj 2 - (q 2 /2mj - qv Fi ) 2 ~\ 
- {q 2 /2rrii + qv F i) 2 \ ' 



IlS(g, W ) = ^ln 
nq 



where the principal value of logarithm (i.e., —n < Im[ln(z)] < tt) should be taken for complex 
frequency. In evaluating U. (q,uj) for real frequency, the usual retarded limit uj — > uj + i0 + 
is implied. In Eq. (0) and Vpi are the effective band mass and the Fermi velocity, 
respectively. In the long wavelength limit (q — > 0), we can expand the irreducible ID 
polarizability as 

II°( g ,a;)«^ 2 f (8) 

where vpi = kpi/rrii is the Fermi velocity with a Fermi momentum k Fi = irrii/2 for the i-th 
ID component with rij being the ID density in the z-th component. Since RPA is known to 
be valid in the long- wavelength limit (q — > 0), the limiting forms for the polarizability would 
be sufficient for our analytic calculations since we are interested in the leading order wave 
vector dependence of the collective modes of the system. [For our numerical results we use 
the full ID II(g, u) as given in Eq. (|7|).] Note that II, v F , k F , etc., refer to the noninteracting 
or the bare ID system, and are therefore perfectly well-defined. 

The expression for the ID polarizability n (g, uj) given in Eq. (|7|) assumes that the sys- 
tem is pure and free from any impurity scattering. In the presence of impurity scattering 
characterized by a level broadening or damping 7, the polarizability Il 7 can be evaluated 
by means of the standard perturbation theory [Tj| to include the impurity scattering in- 



duced vertex correction. The impurity scattering effects are introduced diagrammatically in 
the RPA by including impurity ladder diagrams into the electron-hole bubble consistently 
with self-energy corrections in the electron Green's function. Since the exact expression 
for TLy(q,u) within this diagrammatic approach is complicated, we use in the numerical 
calculation a particle-conserving approximation for arbitrary values of q and uj, given by 



Mermin ||18|| . The Mermin expression is, in fact, equivalent to the diagrammatic result at 



long wavelength as long as the density of states renormalization by impurity scattering is 



negligible. In this relaxation time approximation [18|, with the impurity scattering induced 



level broadening 7 = l/(2r), the polarizability is given by 



n / \ _ (uj + i-f)Il (q,uj + i-f) 
h[q,UJ) u, + vy\no{q,u, + vy)/IIo(q,0)Y 
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For small q and u, the polarizability within the Mermin formula becomes 



n 7 (?,w) 



nq 



muj(u + ij) 



(10) 



This approximate formula gives the same long wavelength and low frequency diffusive be- 
havior as one gets from the diagrammatic approach ||15fl . In this paper, we take the impurity 
scattering induced broadening 7 as a constant phenomenological parameter, which, for ex- 
ample, could be taken from the experimentally measured carrier mobility \i = e 2 r/m, with 
7 = l/(2r). 

The condition for the existence of a collective mode is given by the zeros of |e|, the 
determinant of the dielectric tensor defined by Eq. ([TJ). For a two-component system without 
any tunneling, we have within RPA: 



1 - V(q) \u° u (q,u) + n° 2 (g,u;)l + \V(q) 2 - U(q) 2 } U° u (q, u)U° 22 {q, u) = 



(11) 



with i,j = 1,2 being the two charge components in confined wires 1 and 2. This equation 
indicates that the sign of the charge components forming the plasma has no effect on the 
dispersion relation of the collective modes of a two-component system. 

The analytical formulae for the long wavelength plasmon dispersion can be obtained 
even in the general situation where the two charge components have different Fermi wave 
vectors and Fermi energies (by virtue of having, for example, different effective masses m 12 
and densities rz 12 ) • In this section we study the collective modes in this general situation 
(i.e., d 7^ 0, rii 7^ n 2 , and mi 7^ m 2 ) although we should emphasize that such a completely 
general situation is not physically particularly applicable (because in most existing quantum 
wire systems either one has mi = m 2 and d 7^ or mi 7^ m 2 and d = 0, corresponding 
respectively to a bi-wire system or an EHP system), and in Sec. IV, we investigate, in detail, 
the modes in a photoexcited homogeneous (i.e., d — 0) ID EHP (d — 0, n\ — n 2 , mi 7^ m 2 ) 
where the electrons and holes are not spatially separated, but have different effective masses. 

First, we consider the plasmons of a symmetric two-component system having equal 
densities {n\ = n 2 = n), masses (mi = m 2 = m), and a finite inter- wire separation d, and 
neglecting impurity scattering effects (7 = 0) for simplicity. In this case IT^i = IT 22 and Eq. 
([TT|) becomes 



l-2V + (q)U° u (q,u) x 1 - 2V-(q)U° n (q, u 



0. 



(12) 



where V± = [V(q) ± U(q)]/2. By solving Eq. fljjP we have two branches of collective modes 
corresponding to the in-phase OP mode u> + and the out-of-phase AP mode u;_ (the in-phase 
uj + and out-of-phase ui- nomenclature applies to the situation when both components are 
electrons or both are holes; otherwise, i.e., for an electron-hole system, uj + is the out-of-phase 
mode and a>_ the in-phase mode): 
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A ± (q)[v F q + E(q)] 2 -[v F q-E(q)f 
A±{q) - 1 



(13) 



where E(q) = q 2 /2m and A±(q) = exp[7rg/(2mV±)]. The optical mode ui + (acoustic mode 
uoJ) indicates the symmetric in-phase (antisymmetric out-of-phase) density oscillation of the 
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coupled system assuming the two carriers have the same charge. In the long wavelength limit 
q — > 0, Eq. fljjD becomes 



w±(g -»• 0) = q 



v 2 F + -w F V±(g -> 0) 



7T 



1/2 



(14) 



where the asymptotic forms of the Coulomb potential V±(g) as q — > are given by 

2e 2 

-> 0) = — [K (\qa\) + 1.97269...] , 



2e 2 r 1 

V(q -»• 0) = / rfxln 

e b Jo 



3 

'1 - x) (2 + cos(27rx)) + — sin(27ra) 

2tt 



(15) 



where = [x 2 + (<i/a) 2 ] 1//2 = x + d/a] for the system aligned in the x-z (x-y) plane. 

As q — > 0, V r + (g)diverge as — ln(ga), but the leading order term in V-(q) is finite and depends 
only on the spatial separation d. Therefore, the leading order term in the ID OP dispersion 
u) + (q) is enhanced by a factor y/2 compared to the one-component system (corresponding 
trivially to the change in density from n to 2n) and behaves as q\ ln(ga)!] 1 / 2 as q — > 0. The 
dispersion of the acoustic mode ui-(q) shows a purely linear behavior cu_(g) = v_q for q — ► 
with the velocity t>_ = Vf[1 + AV-(q = 0)/ (tivf)} 1 ^ 2 as g — > 0, and lies always above the 
single particle Landau damping region because V- > vp. The spatial separation of the two 
wires produces a weak second order correction to the OP dispersion in the two-component 
system. For d = 0, the u;_ mode becomes degenerate with the electron- hole continuum (i.e., 
v _ = vp) in the long wavelength limit, and the u + mode becomes that of the corresponding 
one-component system with a total density iV = In. As the spatial separation increases the 
optical (acoustic) plasmon mode frequency decreases (increases) slowly at long wavelengths. 
For the system with an infinite separation, the two modes become uncoupled, and have the 
same dispersion as the ID plasmon in the one-component system with a total density n. 

In Fig. |l] we show our calculated collective modes dispersion for different spatial 
separations between the two identical quantum wires; (a) d = 0.2a^ and 2.0a^ with 
n = 0.6 x 10 6 cm -1 , and (b) d = 2.0a_e and 4.0a B with n = 10 6 cm -1 , where a B (~ 100 A) is 
the Bohr radius of the GaAs system and we use this as our length unit throughout this paper. 
In Fig. |I] we use the parameters: m = 0.067m and the wire width a = a B . As discussed 
above, we find that the plasmon dispersion strongly depends on the interwire Coulomb cor- 
relation determined by d. In these figures we find that for larger spatial separation the 
energy difference of u± modes is small because as d increases the interwire Coulomb corre- 
lation decreases. The important point to note is that both the u± modes remain undamped 
upto very large wave vectors due to the suppression of ID SPE, and therefore experimental 
observation of u± charge density modes in coupled bi-wire systems should be relatively easy 
via light scattering or far-infrared spectroscopy. 

Now we consider the collective mode spectrum in a general system where the two com- 
ponents may have different densities and masses. The general situation with d = is of 
special experimental relevance because it describes the photoexcited ID EHP system - 
we therefore defer the discussion of the d = ID EHP system to section IV. Since there 
is a low energy gap in the ID SPE, we look for undamped plasmons in two regions: the 
high frequency region (u > qvpx, qvF2) and the low-frequency region (qvpi > oj > qvpz)- 



S 



(We choose Vfi > Vf2 throughout this section without any loss of generality.) 
wavelength limit (q — > 0), we solve Eq. ( |TT|) using Eq. (||) to get 



In the long 



uj + (q) = qv Fl 



1 + ^v l(s ^0) 



1/2 



(16) 



r 8 a l 1 / 2 
co_(q) = qv F1 a + ^(<? ^ °) ' ( 17 ) 

where a = Vpz/vpi = {n 2 /ni) (mi/m 2 ) < 1. In the long wavelength limit, these asymptotic 
expressions for uj± are valid in both the high frequency (u > qvpi) and the low frequency 
(qvpi > uj > qvp2) regimes. The condition for the existence of an undamped AP mode 
in the high frequency regime (uj > qvpi) is given by v^(a, d) > vfi, or 

a+ 8 ° VL(g->0)>l. (18) 
7rt>i?i 1 + a 

We find that the condition defined by the inequality [Eq. (|18D 1 leads to the conclusion 
that within RPA can exist as an undamped mode either for a large spatial separation 
(and arbitrary a) or for almost equal Fermi velocities, i.e., a ~ 1, (and arbitrary d). As 
a decreases the system should have a large spatial separation to have an undamped high 
energy AP mode. For a = 1 (e.g., equal mass and density), we find that the acoustic mode 
lies in the high frequency regime and is always undamped for arbitrary d. This is also true 
in 2D bilayer systems with finite layer separations, where the AP mode is undamped as long 
as a = 1 and d ^ 0. The condition for the existence of the AP mode in the low frequency 
regime (qvpi > > qvF 2 ) is given by vfi < v-.(a,d) < Vfi- Due to the energy gap in 
the ID SPE the AP mode is undamped in this regime. The existence of the undamped AP 
mode in the low frequency regime for arbitrary a and d is a unique phenomenon in the ID 
two component system arising entirely from the non-existence of low energy SPE in ID. In 
2D or 3D two-component system the low frequency regime is necessarily accompanied by 
Landau damping due to the presence of the low energy SPE continua. 

In Fig. H we show our numerically calculated longitudinal collective modes (for different d 
and a) with a fixed density n\ = 0.6 x 10 6 cm -1 , effective mass vii\ = 0.067m e (corresponding 
to the effective electron mass of the GaAs), and a wire confinement a = a#. In these 
figures the higher (lower) energy mode oj + {ujJ) denotes the optical (acoustic) mode. In 
Fig. |2](a) we use the parameters: the spatial separation of the two wires d = 0.5a# and 
a = {n,2/ni) (mi/ 7712) = 0.6. (Note that a determines the density of the component two 
for a given n\ and mi, 777.2. For example, if rri2 = mi we have 112 = 0.6ni.) With these 
parameters the velocity of the acoustic mode u_ exceeds the Fermi velocity of the component 
one (f_ > Vpi) and lies well above the SPEs of both components (high frequency regime). 
In Fig. 0(b) we choose the parameters which limit the t>_ between vp\ and Vp2- d = 0.2a# 
and a = 0.2. (Note that this requirement restricts d to unrealistically small values.) Thus, 
the undamped AP mode with these conditions lies in low frequency regime (vp2 < v - < Vfi) 
which corresponds to the regime between two SPEs. The important point is that in both 
cases ("large" d and a, and "small" d and a) the AP is undamped in ID in contrast to 
higher dimensional systems. 
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Next, we investigate the plasmon modes of the two- component system in the presence 
of impurity scattering. We calculate the plasmon dispersion of a spatially separated ID 
two-component system with equal masses and densities (a = 1) in the presence of impurity 
scattering (characterized by a constant level broadening 7) . With the approximation for the 
polarizability function for small q and u defined by Eq. ( |I0| ) the plasmon dispersion is given 
by 

^~-il + ^-t + ^ V± ( q ), (19) 

where u+ (u-) indicates the OP (AP) mode. Note that for 7 — > 0, Eq. (|19D recovers the 
result previously obtained for the pure two-component system. Thus, in the presence of 
impurity scattering the plasmon modes acquire a wave vector dependent imaginary part 
which corresponds to a plasma damping, uj±(q) = uj p ±(q) — ia±{q). The plasmon becomes 
overdamped at small q as it becomes completely imaginary for small q, i.e., the plasmon 
modes are totally overdamped below the critical wave vectors q c ± and exist only for wave 
vectors q larger than these critical values. (Note that the plasmon broadening ct±(q) is, in 
general, determined by the imaginary part of the complex plasmon mode, and is different 
from the level broadening 7.) For q > q c ± the plasma dampings are roughly comparable to 
the half of the impurity scattering, i.e., u±(q) ~ 7/2. From Eq. (|T9|) we find the critical 
wave vectors (below which the plasmon is no longer a well defined collective mode) as 

K 

Qc+ ~ y/2\ ln(ira/\/2)|' 
K 

Qc- = , , 

where K = r y^Jmei ) /(8ne 2 ). As the impurity scattering rate 7 increases the critical wave 
vectors q c ± also increase. The q c+ (just as u + itself) depends weakly on d, but g c _ (just 
as ujJ) is strongly affected by d. As d —>■ the intra- and inter-wire Coulomb interactions 
become the same, V^(q) — ► 0, and consequently g c _ — > 00. Thus, for the system with a 
large impurity scattering and a small spatial separation we have only the OP mode in the 
long wave length limit. In general, we have g c _ > q c+ from Eq. (EH), which means that the 
AP mode is easily affected by impurity scattering (and is therefore more difficult to observe 
experimentally.) 

In Fig. |3] we show the plasmon dispersion calculated numerically by using the full polar- 
izability, Eq. (|pD , in the presence of finite impurity scattering. We find zeros of the complex 
dielectric function of the ID two-component system in the presence of impurity scattering, 
i.e., e(q,u) = 0, to obtain the plasmon dispersion curves. In Fig. |] the curves with lo > 
give the plasma frequency or the real part [w p (g)], and those with to < give the plasma 
damping or the imaginary part (a) of the complex zero solutions [&>(g)]. The figure shows 
that the plasmons are overdamped below critical wave vectors q < q c ±. The critical wave 
vectors, q c ±, below which the plasmon does not exist due to impurity scattering effects, 
depend on the density and the impurity scattering rate. We note that the overdamping of 
the plasmon mode occurring for small q is a direct consequence of the diffusive nature of 
the electron dynamics. We emphasize that this long wavelength impurity scattering-induced 
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(20) 



suppression of collective modes is more severe in ID than in 2D or 3D (in 3D, because the 
OP is gapped, it is hardly affected by impurity scattering.) 

Finally, we study the dynamical structure factor, S(q,u) oc Im[e(g, cu) -1 ], which gives 
a direct measure of the oscillator strength of the density fluctuation spectrum, and many 
experimental probes such as inelastic electron and Raman scattering spectroscopies are 
directly related to the dynamical structure factor ||19|| . Equivalently, the dynamical structure 
factor is a measure of the spectral weight carried by the particular collective mode of the 
electron liquid. An undamped plasmon shows up as a 5-function peak in S(q,ui) indicating 
the existence of a simple zero of e(q,u). When both Re[e] and Im[e] become zero [i.e., 
e(q, to) = 0] the imaginary part of the inverse dielectric function Im[e(g, u;) -1 ] is a 5-function 
with the strength 



where ui±(q) are the plasmon dispersions for the system. The damped plasmon (a ^ 0), 
however, corresponds to a broadened peak in S(q, uj) - for larger broadening, as for q < q c 
in the system with impurity scattering, the plasmon is overdamped and there is no peak 
in S(q, to). At long wavelengths (q — > 0), the optical plasmon u + mostly exhausts the 
/-sum rule and carries almost all the spectral weight. W(q) approaches zero as q — ► 
because the ID plasma frequency vanishes. Note that the small weight of the AP mode 
at long wavelengths makes it particularly susceptible to collisional damping effects in ID 
two component system. For the system with the same effective mass and a large spatial 
separation, the acoustic plasmon mode has an energy comparable to that of the optical 
plasmon mode, which gives the enhanced spectral weight associated with the AP mode. 

In Fig. f| the Im[e(g, uj)' 1 ] is plotted for different wave vectors as a function of the 
frequency for two different systems. In Fig. ^(a) we choose the parameters such that the 
AP mode lies in the low frequency regime: rii = 1.0 x 10 6 cm _1 , a = 2a#, a = 0.2, and 
d = 0.5aB- In Fig. |](b) we choose the parameters such that the AP mode lies in the 
high frequency regime: ri\ = 0.6 x 10 6 cm _1 , a = as, a = 0.9, and d = as- The insets 
in Fig. |4] show the weight W(q) of the plasmon modes. For pure system the plasmon 
peaks become 5-function peaks with weights given in the inset, but with impurity scattering 
rate 71 = 0.1-E/i and 72 = 0.1-E/ 2 , we show broadened peaks giving the damped plasmon 
modes. As shown in Fig. ^ the spectral weight of the acoustic mode is comparable with that 
of the optical mode at finite wave vectors, but not in the long wavelength limit. With the 
enhanced spectral weight the acoustic plasmon mode in ID should be more easily observable 
than in higher dimensions and may significantly affect physical properties of the quasi- ID 
two-component systems. The results presented in this section indicate that the AP mode 
should be observable in double wire structures under experimentally realizable conditions. 

III. BI-WIRE WITH WEAK TUNNELING 

In this section we analytically calculate the collective mode dispersions in bi-wire struc- 
tures (a = 1, i.e., rii = n 2 and m x = m 2 , and d 7^ 0) in the presence of significant interwire 
quantum tunneling, all earlier work in the previous section having dealt with the zero tun- 
neling limit. Since our focus is on understanding tunneling effects in this section, we put 



11 



a = 1 with no loss of generality. In the presence of interwire tunneling, the electron energy 
eigenstates E± should be used [g0| as the basis set rather than the wire index which is no 



longer a good quantum number. The energy levels E± = e(k) ± t, where e(k) = k 2 /2m is 
the parabolic one electron ID kinetic energy in each wire and t is the tunneling strength, are 
the usual symmetric and antisymmetric one electron eigenstates in the presence of tunneling 
with the single particle symmetric-antisymmetric (SAS) gap given by As as = E + — E_ = 2t. 
In the SAS representation the collective mode spectra become decoupled by virtue of the 
symmetric nature of our bi-wire system (i.e., both wires identical with equal electron den- 
sity), and the collective density fluctuation spectra are given by the following two equations 
for the in-phase and the out-phase plasmon modes uj± respectively: 

e+(q,u) = 1 - V+(q) \a++(q,u) + n__(g,w)] = 0, (22) 

and 

e_(<Z,u;) = 1 - V-(q) [Tl + 4q,u) + n_ + (g,w)] = 0, (23) 

where V±(q) = V(q)±U(q) with V(q) and U(q) being respectively the intrawire and interwire 
Coulomb interaction matrix elements given in Eqs. (@) and @. The U a p(q,u}) in Eqs. (|22|) 
and (f2"3"|), with (a, (3) = (+,—), are the noninteracting SAS polarizability functions within 
our RPA theory: 

J (2ir) 2 w + E a (k + q) - ^(k) 

where f a ^ are Fermi occupancy factors. 

Solving Eqs. (p^)-p^) we obtain the collective density fluctuation spectra of the coupled 
bi-wire system. In the absence of tunneling, t = 0, one has E + = E_ = e(k), and one then 
recovers in a straightforward fashion the results we have in Sec. II, i.e., the in-phase optical 
(u+ ~ q\ ln(ga) \ 1 ^ 2 ) and the out-of-phase acoustic (u- ~ q) plasmons of a bi-wire system 
without any electron tunneling. It is, in fact, straightforward to obtain analytically [from 
Eqs. (|22])-(p^)1 the long wavelength (q — > 0) plasma modes of the coupled bi-wire system 
including effects of interwire tunneling. We obtain in the long wavelength limit the following 
results: 



co 2 + (q^0) = q 



+ -(b + + b„)v F V+(q) 



(25) 



ul{q - 0) = A 2 SAS + -^{b + - 6_)VL(g = 0)A 5A5 , (26) 

7T 

where Vp = kpjm is the Fermi velocity and b± — (1 ± t/Ep) 1 ^ 2 for n > n c = AmAsAs/^ 2 
when both symmetric and antisymmetric levels are occupied (i.e., the ID Fermi energy 
Ep > As as), an d b + = 2, 6_ = for n <n c when only the symmetric level is occupied (i.e., 
Ep < A S as)- 

The most important qualitative feature of the plasmon dispersion in the presence of 
interwire tunneling is that even though the in-phase OP mode (uj + ) depends weakly on 
tunneling, the out-of-phase mode (u;_), which is purely acoustic in the absence of tunneling 
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~ q if As as = 0), develops a plasmon gap at q = in the presence of nonzero tunneling. 
The plasmon gap A = u)-(q = 0) depends nontrivially on the ID electron density n, the 
interwire tunneling amplitude t, and the interwire distance d. It is easy to see from Eq. ( |26|) 
that this plasmon gap A, has the following behavior 

A ~ A s as or \J A s as, (27) 

depending on whether the interwire tunneling is strong or weak. It should be emphasized 
that the strikingly non-intuitive A S '2 S dependence of the collective mode gap (on the square 
root of the single particle gap) is purely a Coulomb interaction effect, which dominates 
the collective excitation spectra in the weak tunneling situation. We mention the curious 
phenomenon that interwire tunneling converts the u;_ mode into an effective optical mode 
by opening up a long wavelength plasmon gap (for t = the U- mode is the AP mode). 



IV. PHOTOEXCITED ONE DIMENSIONAL ELECTRON-HOLE PLASMA 

In sections II and III we discussed a double wire system where the two one dimensional 
charged components are separated by a distance d. Now we consider the possibility of acous- 
tic plasmons in a photoexcited homogeneous (i.e., d = 0) ID EHP where the electrons and 
holes are not spatially separated. This is, in fact, the original context in which the AP 



mode was first discussed in the solid state plasma more than forty years ago. The photoex- 
cited EHP in semiconductors is also the system in which most of the early investigations 



lT| of acoustic plasmons were carried out. Our reason for singling out the ID EHP for a 
detailed investigation in this section (in principle, it is just the d = limit of Sec. II) is 
the existence of ID EHP in real samples from several groups where a search for ID uj± 



modes via light scattering spectroscopy should be successful. 



It has been known for a long time [|10j that a two-component EHP allows for the existence 
of two collective modes, one of which (the low frequency one) is linear in wave vector or 
acoustic at long wavelengths, and the other mode (the high frequency or the optical plasmon) 
is essentially the combined collective charge density oscillation of the whole system being 
therefore qualitatively similar to the usual plasmon mode in a single component plasma. 
In an electron-hole two-component plasma, as appropriate in a photoexcited semiconductor 
system, the AP corresponds to the in-phase collective density fluctuation excitation of elec- 
trons and holes, and the OP corresponds to their out-of-phase motion with respect to each 
other. [If the two components are both electrons or both holes, for example, in a bi-wire (as 
considered in Sec. II of the paper) or bilayer or a two-band system, then the AP signifies 
the out-of-phase collective mode and the OP the in-phase motion - in general, the acoustic 
mode is the neutral mode at long wavelengths and the optical mode the charged mode.] In 



spite of compelling theoretical arguments ||i0| , |22| , |23l for the existence of "quantum" AP in 



degenerate two-component homogeneous electron-hole solid state plasmas the experimental 
observation of low temperature AP in two and three dimensional homogeneous electron-hole 
semiconductor plasmas has been very severely hindered by Landau damping effects. It turns 
out that in two and three dimensional degenerate homogeneous EHP the AP is necessarily 
Landau damped (by decaying into single particle excitations) and is in fact unobservable 
as a well-defined mode unless the effective mass ratio between the two species is extremely 
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large |10| , p2| , |23| . This has prevented an unambiguous observation of a quantum AP in 
degenerate solid state EHP although the corresponding classical AP (in the non-degenerate 
high temperature electron-ion plasma) was experimentally observed pll almost forty years 
ago. 

In this section we show that the elusive quantum AP in the homogeneous (i.e., d = 



0) EHP should be relatively easy to observe in the degenerate ID EHP [^TJ confined in 
semiconductor quantum wire structures due to the severe suppression of single particle 
excitations in one dimensional systems. In particular, within the RPA which should be 
quantitatively valid in one dimension by virtue of the smallness of vertex corrections, the 
AP mode in a ID EHP in GaAs quantum wires should be completely undamped upto a rather 
large wave vector q > kp, in contrast to two/three dimensional EHP where the AP is damped 
in GaAs and is essentially unobservable. We believe that resonant inelastic light scattering 
spectroscopy, which has already been successfully used || in observing the ID plasmon in 
single component (n-doped) GaAs quantum wires, is the ideal tool in the search for the 
AP in photoexcited (undoped) two component EHP in GaAs quantum wires. Theoretical 
results presented in this section show that the ID in-phase AP mode should be readily 
experimentally observable in resonant inelastic light scattering experiments performed on 
photoexcited ID EHP [2~I| in narrow GaAs-AlGaAs quantum wires. 



We put d = in the formalism developed in Sec. Ill of the paper and allow for the possi- 
bility of different masses m e and taking equal densities n e = = n of the electrons and 
holes. Below we present and discuss our results for collective modes in ID EHP both in the 
RPA and in the Hubbard approximation (HA) where one includes a local field correction to 
the RPA non-interacting polarizability. The inclusion of local field corrections (HA) beyond 
RPA distinguishes this section from Sec. II where strict RPA was used — our motivation 
for using HA (as well as RPA) is to make our results more realistic for the experimental sys- 
tems where local field corrections may very well be important at the experimentally feasible 
densities and wave vectors. The dynamical dielectric function in the HA is given by 

e E (q,u) = {l- V(q)[l - G e (g)]n°(g,^)} {l - V(q)[l - G h (q))lP h (q,u)} 

-V{qfW e {q,uj)W h {q^ (28) 

where G e< h(q) is a simple static local field correction to the RPA and is given by G e ^(q) = 



V{y<l 2 + k 2 Feh ) /2V(q), where k Fe ,h are the Fermi wave vector of the electron and hole re- 
spectively. The calculation in the HA is similar to that in the RPA with the RPA dielectric 
function [Eq. ([T])] being replaced by the HA dielectric function [Eq. ( p8|) 1. By putting d = 
we have V_(q) = 0, V + (q) = V(q). Thus, in the long wavelength limit q — > 0, we have the 
plasmon mode dispersions within RPA from Eqs, (|i"6|) and ([FT 



u+ (q) = qv Fe 



7lV Fe 



1/2 



(29) 



\q) = y/aqv Fe , (30) 



where a = vph/v Fe corresponds to the mass ratio m e /mh since n e = n h . As the mass ratio 
decreases (a — * 0), the OP mode becomes the plasmon mode of the one component 
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system. Since < a < 1 we have qvph < ^-(q) < QVFe- The AP mode within RPA lies 
always in the undamped region (low frequency regime) between two SPEs. This undamped 
ID AP mode uo_ in the long wavelength limit is fundamentally different from the AP mode 
of the higher dimensional EHP since the presence of the low energy Landau damping region 
in 2D or 3D gives rise to the complete damping of AP mode, unless a is very small — in the 
ID EHP the AP mode is undamped within RPA for arbitrary values of a (a < 1 by choice). 

In Figs. [5] -|7| we show the numerically calculated plasmon modes and spectral weights 
for different densities in the photoexcited GaAs quantum wire ID EHP within both RPA 
and HA. In Figs. ^] -|7] we use the parameters corresponding to GaAs: m e = 0.067m , 
m h = 0.45m , e b = 11, and the confinement width of the wire a = as- We also use a level 
broadening 7^ = Q.lEp e> h due to the impurity scattering in calculating the mode spectral 
weights. Note that in the long wavelength limit the RPA AP modes are undamped (except, 
of course, for impurity damping) for all densities because of the absence of long wavelength 
Landau damping. In general, similar to the corresponding higher dimensional situations 
the local field effects in the HA reduce the plasma frequencies |15| compared with the RPA 
results, in particular, in the low density system. Since the uj + mode exhausts the /-sum 
rule in the long wavelength limit, the local field correction has relatively little effect on the 
uj + mode. The uj_ AP mode is, however, very strongly affected by local field effects. Fig. 
|5] shows the results for the EHP density n e = = 10 6 cm -1 : At this density we find the 
remarkable result that the co>_ AP mode is completely suppressed by local field effects as 
it gets pushed into the SPE Landau damping continua. Thus, in Fig. |5](a) the AP mode 
dispersion within HA does not show up because u- lies completely in the Landau damping 
continuum. Fig. 0(b) shows that within RPA (solid lines) the cj_ mode (low energy peaks) 
carries an appreciable spectral weight which is comparable to that carried by the OP mode 
uj + (high energy peaks), however, within HA (dashed lines) the uo + mode actually carries all 
the spectral weight (we find only the high energy uj + peaks in the HA with the low energy uj- 
peaks being completely suppressed.) As density increases the (completely suppressed) oj- 
mode in HA reappears in the low energy regime as the AP mode comes out of the Landau 
continuum. In Figs. ^ and [7| we show the results for density n e = n/j = 1.2 x 10 6 cm -1 
and n e = rih = 1.5 x 10 6 cm -1 , respectively. At the density n e = = 1.2 x 10 6 cm -1 
the o;_ mode reappears below the SPE^ in the HA. This oj_ mode lying below SPE^ is 
undamped due to the low energy SPE gap in ID and is a characteristic feature of ID EHP. 
The spectral weight of the undamped uj- mode in the low energy region (at higher density) 
is comparable to that of the uj + mode even within HA [Fig. §](b)]. At very high density 
(Fig. D n e = % = 1.5 X 10 6 cm -1 the local field corrections are weak, and we find that 
the cu_ mode (in the HA) lies between SPE e and SPE^ and is undamped in the long wave 
length limit. The calculated spectral weight in the HA is comparable to that in the RPA in 
these high density systems. One conclusion of our HA-based calculations is that one should 
use higher EHP densities for an unambiguous observation of undamped ID AP — this is 
consistent with our RPA results also, where Landau damping is absent, because in RPA 
higher carrier densities imply quantitatively weaker effects of impurity level broadening. 

Given the unique nature of the SPE continua (and the strong low energy suppression 
of Landau damping) in the ID system, the question naturally arises about whether the 
collective mode behavior in ID EHP is fundamentally different from that in 2D and 3D 
systems, where the AP mode is usually unobservable due to its unavoidable decay via Landau 
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damping to the electron(i.e., the lighter mass component) SPE unless a is very small (i.e., 
the mass ratio very small, m e <C m^) which cause the Landau damping to be small, but still 
non-zero. In the ID EHP, however, there are three (as opposed to only one in 2D and 3D 
EHP) Landau damping-free regimes in the uj-q space — these are [see Fig. 5(a) or 6(a) or 
7(a)] the high energy-low wave vector regime above the lighter mass SPE continuum (which is 
present in 2D and 3D systems also), the low energy regime below both the SPE continua, and 
the intermediate energy — low to intermediate wave vector regime in the SPE gap between 
the SPE continua of the two components (we emphasize that Landau damping of collective 
modes is allowed in 2D and 3D systems in the last two regimes even in the RPA because 
there is no low energy gap in the SPE continuum in higher dimensions). To emphasize the 
qualitative difference between ID EHP and its higher dimensional counterparts with respect 
to the Landau damping characteristics of the collective modes, we provide in Figs. 8 and 9 
our calculated RPA results for collective mode dispersion and damping in 2D and 3D EHP 
respectively, choosing the same effective mass, lattice dielectric constant, level broadening, 
and "approximately equivalent" carrier densities. One can see that in Figs. 8(a), the 2D 
EHP case, and 9(a), the 3D EHP case, the OP mode is above the SPE continua of both 
components (exactly the same as in Figs. 5-7 for the OP in ID EHP) and is therefore 
not Landau damped whereas the AP mode in both case lies within the electron (i.e., the 
lighter mass) SPE continuum (but above the hole SPE) and is therefore Landau damped 
by decaying into the quasiparticles of the electron component in the EHP. The contrast 
between the RPA results for the 2D/3D AP mode (Figs. 8 and 9) and the ID AP mode 
(Figs. 5-7) is striking in terms of Landau damping — the 2D/3D AP mode is invariably 
Landau damped by being inside the SPE continuum of the lighter carrier component (i.e., the 
higher Fermi velocity, vp e , in our case) where the ID AP mode (within RPA) is undamped 
by being in the intermediate regime in between (the SPE "gap" regime) the SPE continua 
of the two components and is therefore not Landau damped by either plasma component. 
This non-existence of AP Landau damping in the SPE "gap" regime in a ID EHP is a 
fundamental one dimensional effect arising from the suppression of long wavelength single 
particle excitations peculiar to ID. We should mention, however, that, as is obvious from 
Fig. 5-7, inclusion of local field corrections within the HA may push the AP mode down 
into the SPE continua of the higher mass plasma component, i.e., the hole SPE, making the 
AP mode strongly Landau damped even in the ID EHP — our comparison between ID and 
2D/3D EHP collective modes and their qualitative difference are based on the RPA, going 
beyond RPA may very well render the ID AP unobservably Landau damped as happens, 
for example, within the HA in our calculation for n < 1.2 x 10 6 cm^ 1 in Figs. 5 and 6. 

It is very anticlimactic, however, that even within the RPA the qualitative difference 
in the Landau damping behavior of the AP in the ID EHP compared with its higher di- 
mensional counterparts is not dramatically reflected in the mode spectral weight or the 
dynamical structure factor, as can be seen by comparing the RPA spectral weights in Fig. 
5(b), 6(b), 7(b) in the ID EHP with those in the 2D/3D EHP given in Figs. 8(b)/9(b). 
(We emphasize that we use equivalent impurity level broadening, 7 = 0.1E F , in all our 
spectral weight calculations.) In general, the spectral weight carries by the AP mode (the 
lower frequency peak in the dynamical structure factor Im[e _1 ] for each q) is comparable, 
albeit slightly larger, in the ID EHP as that in the 2D/3D EHP. In particular, the AP 
spectral weights in the ID and the 2D EHP seem to be quantitatively similar whereas the 
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corresponding 3D AP mode carries considerably less spectral weight (which is understand- 
able based on the fact that in 3D the OP has a finite gap at long wavelengths in contrast 
with the 1D/2D OP which vanishes at long wavelengths, and therefore in a 3D EHP the 
long wavelength OP mode carries considerably more relative spectral weight than its lower 
dimensional counterparts). The fact that the ID AP carries small spectral weight in spite of 
its lack of Landau damping can be easily understood on the basis of the /-sum rule, which 
is essentially completely exhausted by the OP in the long wavelength limit, leaving rather 
little spectral weight to be carried by the ID AP in the long wavelength. At large wave 
vectors (i.e., away from the long wavelength limit), on the other hand, the AP is Landau 
damped (because it enters the SPE continua) and all collective modes start carrying little 
spectral weights. 

We therefore conclude this section on ID EHP collective modes with two relatively mod- 
est conclusions: (1) In spite of the non-existence of AP Landau damping in the ID EHP 
within RPA, the calculated AP spectral weight is not much larger than the correspond- 
ing 2D situation (it is, however, substantially larger than the corresponding 3D situation), 
and therefore the observability of ID EHP AP should be only slightly more favorable than 
the corresponding 2D case; (2) inclusion of local field corrections in general severely sup- 
presses the AP spectral weight in the ID EHP, and at least within the HA the AP becomes 
unobservable at lower densities. 



V. TWO-COMPONENT LUTTINGER LIQUID 



In this section we briefly discuss, for the sake of completeness, the AP mode of the 
two-component Luttinger liquid, which has two different Fermi velocities vp e and vph corre- 
sponding to the electron and the hole respectively. We make the standard Luttinger model 
approximation of linearizing the single particle energy. This system can be described by the 



2-component Luttinger model [^5] Hamiltonian 



H 



V Fi k 

i,k,a 



a i,k,a a i,k,cr 



i,j,q 



(31) 



where i,j=e,h, and (h,k,a) is the destruction operator of the right- (left-) moving i-th 
component with momentum k and spin a, Vij(q) the Coulomb interaction between i-th and 
j-th component (here Vij = V ee = Vhh = —V e h, and Pi(q) = Pi, a (q) + Pi,b(q) the total density 
operator of the i-th component. 

By standard bosonization methods p9|j25[] we can diagonalize Eq. ( PI|) and find two 
eigenmodes corresponding to the charge density excitations 



<4(?) 



± 



<fv 2 Fe 



1 + a 2 + (1 + a)V(q) 
1 + a 2 + (l + a)V(q) 



4a 



a+ (l + a)V(q) 



1/2 



(32) 



where a = VFh/vFe and V(q) = 2V ee (q)/(7iVF e )] note that V et 
In the long wavelength limit, q — > 0, Eq. fl32|) becomes 



hh 



-V eh in the ID EHP. 
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L0 A 



(q) = v Fe q 1 + (1 + a)V(q) 



1/2 



(33) 



tO-(q) = Vaqv Fe . (34) 

These long wavelength charge density modes defined by Eqs. ( |3"3"D and ( |3"1"D are exactly the 
same as the long wavelength collective modes in the ID EHP within RPA [Eqs. ( p9| ) and 
(|50D1. This can be explained by the fact that vertex corrections to the irreducible polariz- 
ability of the ID system vanish || in the Luttinger model, and therefore RPA and Luttinger 
model results become identical in the long wavelength limit. Even though the charge density 
mode dispersions of the Luttinger model are identical to the RPA collective mode predic- 
tions based on the Fermi liquid theory, the damping properties of the u>- AP mode in the 
Luttinger model is strikingly different form the RPA result. In the previous section we 
found that within RPA the AP mode in the ID EHP exists only at long wavelengths as the 
AP mode gets severely damped by Landau damping at finite wave vectors where it enters 
the electron single particle continua and is Landau damped. However, since there are no 
equivalent SPE Landau damping mechanisms in the Luttinger model the AP mode in the 
Luttinger model seems not to be damped out even at the finite wave vectors, and seems to 
exist at arbitrary wave vectors. 

The above-discussed difference between the RPA (Fermi liquid) result and the Luttinger 
model (non-Fermi liquid) result seems to suggest a direct experimental approach to observe 
the predicted non-Fermi liquid-like behavior of a one dimensional system by studying its 
plasmon damping behavior which, according to the RPA (Luttinger) model, should (should 
not) manifest finite wave vector Landau damping. This is, however, misleading because 
the strict non-existence of any damping of the collective modes at any wave vectors is 
purely a result of the Luttinger model (rather than being a generic Luttinger liquid result), 
arising entirely from the linearization of the single particle energy [the first term in Eq. 
(j31~l)1 in the Hamiltonian. In the generic Luttinger liquid behavior, where one does not 
linearize the single particle energy, higher order processes (which asymptotically vanish 
at low energies and long wavelengths) produce collective mode damping arising from the 
(multiboson) boson-boson collision processes within the bosonized theoretical description. 
Such multiboson collision processes, which are akin to phonon-phonon anharmonic scattering 
processes producing phonon decay beyond the harmonic crystal theory, are irrelevant in the 
renormalization group sense at long wavelengths, but do produce damping of the collective 
modes at finite wave vectors. In this sense, therefore, no fundamental distinction exists 
between the two theoretical approaches since both predict damping of the collective modes 
at finite wave vectors, but not at zero wave vectors. One expects the AP damping in 
the Luttinger liquid theory to be a smooth function of wave vector q, vanishing only in 
the long wavelength limit q — > 0. This seems to be drastically different from our RPA 
damping result where the AP is undamped upto a wave vector q c when it becomes damped 
as it enters the SPE continuum. Going beyond RPA, however, one finds that the AP is 
damped even outside the SPE Landau damping continua due to multipair (i.e., multiple 
electron-hole excitations) productions, and therefore there may not be much qualitative 
difference between RPA and Luttinger liquid theories even for the AP damping properties. 
This issue, however, merits further investigations, and detailed theoretical calculations of 
the AP damping properties in the two-component Luttinger liquid may very well point to 
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some significant differences with our RPA results, which then may lead to a direct method 
of differentiating between Luttinger liquid and RPA results, based on the experimentally 
measured damping properties. Very recently, such a calculation [55| has been carried out for 
a one component clean Luttinger liquid using the short range interaction which, however, 
has to be generalized to the realistic long-range Coulomb interaction in ID quantum wires 
before any firm conclusion on the damping behavior contrast between RPA and Luttinger 
liquid may be reached. We emphasize that the AP damping property being manifestly a 
non- universal property (determined by irrelevant operators such as the band curvature), its 
actual quantitative calculation in the Luttinger liquid theory is problematic. 



VI. SUMMARY 

We have obtained the dispersion and spectral weight of collective charge density exci- 
tations in zero temperature two-component one dimensional quantum plasma confined in 
semiconductor quantum wire structures. In some sense this paper is a two-component gen- 
eralization of our earlier work [TJ| where we studied the zero temperature collective charge 



density excitations of the one component ID system in some details. In a separate earlier 
publication [jnj] we consider a specific example of collective modes in a two-component ID 
electron system (namely, a two-subband ID system) in the context of obtaining quantita- 
tive agreement with the experimental results of ref. |§. In the current paper we refrain 
from repeating any results from these earlier publications |H|,I5] of ours, and the current 



paper, along with refs. |L3| and [ , forms a reasonably complete quantitative description of 



possible low lying collective charge density excitations in ID semiconductor quantum wire 
structures. 

Our main results in this paper are the following: (1) there are two possible collective 
charge density excitations in a two-component ID quantum plasma associated with the 
relative phase fluctuations in the two charge densities — the OP mode is the in-phase (out- 
of-phase) mode when the two components are the same (opposite) kind (both electrons or 
both holes versus one component electron and the other component hole) and the AP is the 
out-of-phase (in-phase) mode in the same situation (in general, the AP is essentially a neutral 
excitation at long wavelengths and the OP is the collective oscillation of the total charge 
density); (2) in contrast to higher dimensional systems where the AP mode is often Landau 
damped because it lies in the SPE continuum of the faster moving charge component, the 
ID AP mode is invariably undamped at long wavelengths (within RPA) due to the severe 
suppression of long wavelength SPE continua in ID; (3) this RPA result of the non-existence 
of long wavelength Landau damping in the ID AP mode is drastically affected by local field 
corrections, which at low density may overdamp the AP mode by pushing it inside the 
SPE continua of both components — the OP mode on the other hand is relatively robust 
with respect to local field corrections since its long wavelength dispersion is fixed by the 
/-sum rule; (4) the AP mode in general carries little spectral weight at long wavelengths, 
and even at finite wave vectors its spectral weight is only slightly enhanced with respect to 
the corresponding 2D two-component situation, leading to the somewhat disappointing and 
unexpected conclusion that in spite of the complete suppression of the long wavelength SPE 
continua in ID systems the neutral AP mode is not substantially easier to observe in ID 
than in 2D/3D systems; (5) finite inter- wire separation in general enhances the AP spectral 
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weight increasing its observability; (6) finite impurity scattering induced level broadening 
gives rise to a critical wave vector below which the collective modes are always overdamped, 
and thus the concept of a long wavelength collective mode does not strictly apply; (7) finite 
inter-wire tunneling produces a long wavelength plasma gap in the AP mode; (8) Luttinger 
liquid theory and Fermi liquid-based RPA theory produce formally identical OP and AP 
collective mode behavior in ID two component systems making it impossible to observe 
any characteristic ID Luttinger liquid signature in the collective charge density excitation 
spectral in one dimension; (9) the AP mode being strongly (and qualitatively) affected by the 
local field corrections (with the OP mode being relatively unaffected), in principle, it should 
be possible to investigate many body effects in ID systems by studying the dispersion and 
damping of the AP mode (the OP mode on the other hand should be relatively insensitive 
to these many body corrections except at rather large wave vectors). 

In our calculations we have concentrated on two distinct types of experimentally realiz- 
able two-component quasi-lD quantum plasmas confined in semiconductor (GaAs) quantum 
wire nanostructures: spatially separated doped quantum double-wire structures and pho- 
toexcited homogeneous ID EHP in a single wire. Both systems are potentially interesting 
from the perspective of ID collective excitation spectra and provide somewhat complemen- 
tary information. Our hope is that our detailed quantitative (both numerical and analytical) 
investigation of the collective mode dispersion and spectral weight in two-component ID 
system will motivate experimental work (involving inelastic light scattering and far infrared 
frequency domain spectroscopies) on the subject searching for the ID AP mode in quantum 
wire structures. 
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FIGURES 



FIG. 1. Calculated RPA collective mode dispersion for different spatial separations between 
the two identical quantum wires; (a) n\ = ri2 = 0.6 x 10 6 cm -1 and (b) n\ = 712 = 10 6 cm -1 . Here, 
the higher (lower) frequency mode corresponds to the OP (AP) mode in each case. Shaded region 
indicates the single particle excitation (SPE) continuum. 

FIG. 2. Calculated RPA collective mode dispersion (a) for the spatial separation d = 0.5a# 
and a = vf2/ v fi = 0.6; and (b) for d = 0.2a# and a = 0.2. The lo + (^-) lines denote the OP 
(AP) mode and SPEi (SPE2) denote the single particle excitation of the component 1 (2). 

FIG. 3. RPA plasmon dispersions of the ID two-component system are shown for various 
impurity scattering rates 7 as given in the figure. Here lo > gives the dispersions of the mode 
and lo < the dampings of the mode. Thick (thin) lines correspond to the OP (AP) mode. 
The overdampings of the plasmons for q < q c ± are clearly seen. We use the parameters: a = 1 
(rii = ?i2 = 0.6 x 10 6 cm -1 and m\ = m,2 = 0.067m e ) a = as, and d = 0.5a#. 

FIG. 4. The RPA dynamical structure factor Im[e(g,w) _1 ] in ID two-component systems as a 
function of frequency u for different wave vectors q. We use the parameters: (a) n\ = 10 6 cm -1 , 
a = 2aB, ot = 0.2, and d = 0, (b) for n\ = 0.6 x 10 6 cm~ 1 , a = as, a = 0.9, and d = as- We use 
the impurity scattering rate 71 = Q.lEpi and 72 = Q.\Ep2 in both figures. The inset shows the 
weights W(q) of the plasmon modes in the absence of the impurity scattering. 

FIG. 5. The calculated (a) plasmon dispersions and (b) spectral weights for densities 
n e = n h = 1-0 x 10 6 cm -1 in the photoexcited GaAs quantum wire EHP within both RPA (solid 
lines) and HA (dashed lines). The SPE^/, denote the single particle excitation of the electron and 
the hole respectively. At this density the AP mode w_ is completely suppressed by the local field 
effects in the HA (no dashed line corresponding to the AP mode exists within HA). 

FIG. 6. The same as Fig. |5| for densities n e = = 1.2 x 10 6 cm -1 . Note that the AP mode 
lo- lies below the SPE/j due to the local field effects. 

FIG. 7. The same as Fig. |5| for densities n e = = 1.5 x 10 6 cm -1 . Note that the AP mode 
lo- lies between SPE e and SPE/j. 

FIG. 8. The calculated (a) plasmon dispersions and (b) spectral weights for densities 
n e = n h = 0.6 x 10 12 cm -2 in the photoexcited GaAs 2D EHP within RPA. 

FIG. 9. The calculated (a) plasmon dispersions and (b) spectral weights for densities 
n e = rih = 0.465 x 10 18 cm -3 in the photoexcited bulk GaAs EHP within RPA. 
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